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We consider hard-core bosons on the kagome lattice in the presence of short range repulsive interactions and 
focus particularly on the filling factor 1/3. In the strongly interacting limit, the low energy excitations can be 
described by the quantum fully packed loop coverings on the triangular lattice. Using a combination of tensor- 
product state based methods and exact diagonalization techniques, we show that the system has an extended Z 2 
topological liquid phase as well as a lattice nematic phase. The latter breaks lattice rotational symmetry. By 
tuning appropriate parameters in the model, we study the quantum phase transition between the topological and 
the symmetry broken phases. We construct the critical theory for this transition using a mapping to an Ising 
gauge theory that predicts the transition to belong to the 0(3) universality class. 

PACS numbers: 


I. INTRODUCTION 

Strong correlations in quantum many-body systems can 
yield novel phases of matter at low temperatures. The 
most prominent example is the fractional quantum Hall ef¬ 
fect (FQHE)' which is characterized by its fractionalized 
excitations^ and topological order^“* that manifests the long 
range quantum entanglement present in the underlying many- 
body ground state. While the FQHE requires strong magnetic 
fields, fractionalized states can occur more generally in cor¬ 
related systems that preserve time-reversal symmetry. Eol- 
lowing the seminal works by Anderson^“^ we know that a 
class of two and three dimensional frustrated spin systems 
can realize paramagnetic ground states, dubbed quantum spin- 
liquids (QSL), that support deconfined fractionalized S = 1/2 
(spinon) excitations as well as various forms of topological 
orders. 

Over the past decades, there has been a very active search 
for model systems stabilizing different types of topologically 
ordered states.A fertile ground to realize such phases 
are frustrated systems in which the geometry of the lattice 
and/or competing interactions prohibit a simultaneous min¬ 
imization of all the inter-particle interactions."’'^ This can 
lead to the suppression of conventional (spontaneous sym¬ 
metry breaking) orders and favor more exotic ground states. 
An associated question of much interest is about quantum 
phase transitions between topologically ordered and conven¬ 
tional symmetry breaking phases. In this regard, systems of 
hard-core bosons with mutual short range repulsions on var¬ 
ious frustrated lattices have attracted much attention because 
they present a set of rich phase diagrams. A particularly 
interesting case is the kagome lattice where, for suitable fill¬ 
ings, a variety of numerical and analytical evidences now 
point to the existence of a topologically ordered Z 2 liquid 
phase over an extended parameter regime.'®’'^ The low en¬ 
ergy physics of such models in the strongly correlated limit 
is generically described by quantum dimer models (QDM)*^ 
on triangular lattices. This latter class of Hamiltonians is 
known to harbor points in their parameter space, the so called 
Rokhsar-Kivelson (RK) points,where Z 2 topological or¬ 
der is present and strong numerical evidence for its stabil¬ 


ity exists.'® The intimate connection, on the other hand, be¬ 
tween such hard-core boson models and S' = 1/2 models 
with XXZ anisotropy*^ also makes these models relevant to 
research in quantum magnetism. With this insight, perhaps 
it is not surprising that there are proposals that the isotropic 
nearest neighbor Heisenberg antiferromagnet can potentially 
realize a Z 2 topologically ordered ground state. 

In this paper, we explore the strong coupling physics of a 
hard-core boson model on kagome lattice with short range re¬ 
pulsive interactions particularly for the filling factor / = 1/3. 
We show that the effective low energy theory is given by a 
QDM Hamiltonian on a triangular lattice with two dimers em¬ 
anating from each site of the triangular lattice. This effective 
theory is thus equivalent to a quantum fully packed loop (FPL) 
modeP^ on a triangular lattice. Such kind of models in the 
classical version have been extensively studied on some of the 
bipartite lattices.The quantum FPL model on a square 
lattice has been studied by Shannon et. al. where they showed 
a correspondence to the XXZ model on the checkerboard lat¬ 
tice in certain easy-axis limits.^® Using a combination of nu¬ 
merical techniques (tensor product states formalism and ex¬ 
act diagonalization on clusters), we analyze this model on the 
triangular lattice. Our numerical analysis strongly advocates 
for a rich phase diagram of the quantum FPL on the trian¬ 
gular lattice consisting of an extended Z 2 topological liquid 
phase as well as a crystalline phase, known as lattice nematic 
(LN),^" that breaks the three-fold rotational symmetry of the 
lattice (but not the translation symmetry). Taking clue from 
our numerical results, we then construct a critical theory for 
a continuous phase transition between the two phases. Un¬ 
like the usual theories of phase transition, this critical theory 
is not written in terms of the low energy long wavelength fluc¬ 
tuations of the LN order-parameter, but naturally in terms of 
“fractionalized” Ising degrees of freedom sitting at the centers 
of the triangles of the kagome lattice. Mapping the problem 
to the language of Ising gauge theory,^' we can isolate these 
critical degrees of freedom- the so called visons^^ (Ising mag¬ 
netic flux-"*), whose condensation then describes the transition 
from the topological liquid to the LN. The order parameter is 
a bilinear in terms of the visons and hence, the above transi¬ 
tion consists of an example of quantum criticality beyond the 
conventional Landau-Ginzburg-Wilson paradigm."''' Our cal- 
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culation predicts that the critical theory belongs to the 0(3) 
universality class. 

The remainder of the paper is organized as follows. We 
first introduce the model Hamiltonian and derive the effective 
model in Section II. We then present the numerical results and 
conclude about the phase diagram in Section HI. In Section 
IV we derive an effective gauge theory description and dis¬ 
cuss the nature of the transition between the liquid and the LN 
phase. We conclude with a brief summary in Section V. 


II. MODEL HAMILTONIAN 


We start by considering an extended Hubbard model of 
hard-core bosons on the sites of a kagome lattice given by 
the Hamiltonian 


n = nt + nv, (i) 

where 

= ( 2 ) 

{hi) 


describes the nearest-neighbor hopping, with amplitude t, for 
the hard-core bosons that are created (annihilated) by b] (b^) 
on the sites of the kagome lattice and 

Tiv = Vi ^ TiiTij -I V2 ^ ruTij + V3 ^ riiTij 

<*J') {{{iJ))) 

(3) 

i 

denotes respectively the first (Vi), second (V 2 ) and third (V 3 ) 
neighbor repulsive interactions among the bosons (n, = b^) 

along with a chemical potential jl that fixes the particle num¬ 
ber. Interesting physics emerges at rational fractional fillings 
p/q (with p and q being mutually prime and p < q). We shall 
restrict ourselves to the specific fractional fillings of bosons 
which are 1/6, 1/3 and 1/2. 

At integer filling 0 or 1 there are only trivial product states 
possible (due to the hard-core constraint). At fractional fill¬ 
ing factors and the presence of longer-range interactions much 
rich phase diagrams emerge. For example, at 1/3 filling 
and strong nearest-neighbor repulsion (Vi), a plaquette or¬ 
dered state'^ is formed, while for 1/2 filling a uniform su¬ 
perfluid persists for all values of V/t. When turning on fur¬ 
ther neighbor interactions given by V 2 and V 3 , additional in¬ 
sulating lobes emerge at different fractional fillings.Many 
of these bosonic insulators at fractional fillings can host inter¬ 
esting quantum phases with or without spontaneously broken 
symmetry. In the following, we focus on the strong coupling 
phases occurring in the 1/2, 1/3 and 1/6 lobes, with a partic¬ 
ular focus on the 1/3 lobe and show that topological as well 
as long ranged ordered phases can emerge. 

To this end we look at the strong coupling limit of the above 
Hamiltonian. For f = 0, and Vi = V 2 = V 3 = 2V, the 



FIG. 1: (Color Online) (a) Mapping from the 1/3 filled bosonic prob¬ 
lem on kagome to FPL on triangular lattice. Each particle on the 
kagome lattice is mapped to a loop segment on the triangular lattice 
that is obtained by connecting the centers of the kagome hexagons. 
The strength of the 1st, 2nd and 3rd repulsive interactions are denoted 
as Vi, V 2 and V 3 respectively. The limit of Vi = V 2 = V 3 = 2V is 
the focus of this paper. Figure (b) demonstrates the allowed lowest 
order processes in f/V (see the text), (c) Shown is the sublattice used 
for the gauge transformation that changes the sign of g in Eq. (5) as 
described in the text. The sublattice is constructed by the b and c 
sites of the shaded unit cells. 


interaction term can be written as 


'Hv = V ^ ^no 


jf! 


{O} 


4VJ 


161/2 


(4) 


where p = p, + 2V is the effective chemical potential and 
no is the number of particles in each of the hexagons of 
the kagome lattice. It is clear from Eq. (4) that for p = 
4V, 81/, 12V, Hy is minimized by having 1, 2, 3 bosons per 
hexagon respectively (or alternatively filling fraction of / = 
1/6,1/3 and 1/2 respectively). Clearly there are many differ¬ 
ent configurations of bosons that satisfy the above constraint, 
however, we note that since the hexagons share sites, the con¬ 
figurations for different hexagons are not completely indepen¬ 
dent. 

An insight to the number of states spanning the ground 
state sector of Hv for the above commensurate fillings can 
be obtained from the one-to-one correspondence between the 
ground state configurations of the bosons and the hard-core 
dimer coverings on the triangular lattice obtained by joining 
the centers of the hexagons of the kagome lattice as shown in 
Fig. 1(a). Thus each site of the kagome lattice lies on a bond 
of the triangular lattice and the presence (absence) of a boson 
on that site can then be identified uniquely with the presence 
(absence) of a dimer on the corresponding bond of the tri¬ 
angular lattice. This immediately shows that at 1/6 filling, 
the number of ground state bosonic configurations allowed 
hy Hv is equivalent to the number of hard-core dimer cov- 
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erings on the triangular lattice which is known to be exten¬ 
sive (~ in the system size, N (an estimate^® based on 

Pauling’s approximation^^ gives ^ 1.34^). Similarly 1/2 fill¬ 
ing can be cast into a 3-dimer (three non overlapping dimers 
emerging from each site) problem m on the triangular lattice'^ 
where Pauling estimate suggests that the extensive degeneracy 
of the ground state is ^ 1.84^. 

In case of 1/3 filling, which is the specific interest of 
the present paper, we obtain the equivalent fully-packed 
loop (FPL) model on the triangular lattice with two non¬ 
overlapping dimers emanating from each site and the dimers 
form non-intersecting loops as shown in Fig. 1(a). Here the 
Pauling estimate shows that the number of loop configura¬ 
tions scales with system size as ^ 1.7^. Thus, in all the 
above cases, for f = 0 , as expected, the ground state is macro- 
scopically degenerate and has a finite zero temperature en¬ 
tropy. Throughout the rest of this work, we shall exploit the 
above equivalence between the bosons and the dimers and 
shall mostly use the language of the dimers, translating back 
to the bosons whenever applicable. 


A. Effective model in the strong coupling limit 

Small but non-zero hopping (t) induces quantum fluctua¬ 
tions that (e.g. in the form of local ring exchange around small 
plaquettes^^"^') can lift the extensive ground state degeneracy 
by quenching the entropy of the classical model {t = 0 ) ei¬ 
ther by spontaneously breaking one or more symmetries of the 
system (quantum order by disorder'^^'^^) or more interestingly, 
by generating a long ranged quantum entangled state that does 
not break any symmetry of the Hamiltonian (quantum disor¬ 
der by disorder'^^’'^^) but can have non-trivial “topological or¬ 
der”. In this work, we shall show examples of both the routes 
taken by the bosons on the kagome lattice. 

To derive the effective Hamiltonian in the strong coupling 
limit, we treat t/V as a small perturbation parameter to obtain 
(to the leading order in f/L)*^ 

Heff = -5^(6“t5/^6“5“+H.c), (5) 


where g = jV and p,q,l and m are the corner sites of the 
bow-tie labeled by a referring to Fig. 1(b). The sum includes 
all bow-ties that are related by the C 3 symmetry. The center of 
the bow tie may or may not be occupied by a boson in case of 
1/3 and 1/2 filling. Furthermore the sign of g in Eq. (5) can be 
altered by using a simple gauge transformation that multiplies 
all configurations with the factor (—where is the 
number of particles in the sublattice shown in Fig. 1(c). 


B. Dimer representation of the effective Hamiltonian 

We now recast Eq. (5) using the particle to dimer mapping 
mentioned in the previous section. In terms of the dimers, the 
effective Hamiltonian corresponds to the kinetic term of the 



FIG. 2: (Color online) The triangular lattice is put on a torus by 
setting periodic boundary conditions along two independent direc¬ 
tions given by the lattice vectors ai and a 2 . Ao denotes the inter¬ 
section point (a link on the lattice) of the two non-contractible loops 
on the torus. There are four topological sectors characterized by the 
doublets (0, 0), (0,1), (1, 0) and (1, 1)P'^ The local Hamiltonian in 
Eq. (7) can not mix configurations from different topological sectors 
hence, block diagonal in the full configuration basis. 


QDM.'’’ 


OL 

where the sum is over all rhombus shaped plaquettes (a) on 
the triangular lattice. 

In addition to the kinetic term (above) the generic QDM 
also includes a potential term (known as Rokhsar-Kivelson 
(RK) potential),'* which energetically favors configurations 
with plaquettes having parallel dimers. Such potential terms 
are representative of higher order (four boson) terms that can 
be generated in the strong coupling expansion in t/V of the 
underlying boson model. However, here we simply use this 
term as a free tuning parameter in the model. The full Hamil¬ 
tonian of the QDM is written as 

Hrk = -5E(|//)(^ 

O' 

Ol 

where positive (negative) I^k denotes repulsive (attractive) 
interaction between the parallel dimers in a given plaquette. 

This generic form of the above Hamiltonian was first pro¬ 
posed on a square lattice by Rokhsar and Kivelson'* (RK) in 
context of high temperature superconductivity. At a special 
point when g — Vrk^ known as RK point, the spectrum con¬ 
tains zero energy ground state with the wave function given 
by the equal weight superposition of all allowed dimer con¬ 
figurations. The RK point features a U{1) resonating valence 
bond (RVB) liquid phase with algebraic decay of dimer cor¬ 
relations. 

However, on non-bipartite lattices, the scenario changes 
drastically.'®’"'® The RVB phase in this case fosters a gapped 
Z 2 dimer liquid in an extended parameter regime (including 
the RK point) with exponentially decaying dimer correlations 
and is characterized by a topological order in the form of four¬ 
fold ground state degeneracy for a system on a torus (as de¬ 
scribed in Eig. 2). The solid phases, on the other hand, spon¬ 
taneously break various lattice symmetries of the Hamiltonian 
and thus have long range dimer-dimer correlations. 


H.c. 


,(7) 


+ H.c.) , ( 6 ) 
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FIG. 3; (Color online) The phase diagram of the dimer models at 
different fractional fillings /. The important details about each of 
them are given in the text. The focus of the present paper is at / = 
1/3. 


In context of the present work, as noted in the last sec¬ 
tion, we mention that the low energy effective Hamiltonian in 
the strong coupling limit generally assumes the RK form irre¬ 
spective of the filling fraction of the bosons (or the equivalent 
dimer models: QDM, FPL or the 3-dimer model respectively 
for / = 1/6,1/3 and 1/2) considered, albeit with important 
implications for the stability and nature of both the liquid and 
the solid phases which is summarized in Fig. 3. With this gen¬ 
eral formulation we now specialize to the physics of the Mott 
lobe for 1/3 filling. 


III. NUMERICAL CALCULATIONS AT 1 /3 FILLING 

We use a combination of tensor product states (TPS) for¬ 
malism and exact diagonalization (ED) methods to obtain the 
phase diagram (Fig. 3 middle line) of the FPL model on the 
triangular lattice. Necessary details about implementing the 
numerical methods for different clusters are furnished in the 
following subsections facilitating a systematic analysis of our 
model. We take the value of p to be 1 which is a convenient 
choice for further numerical calculations. 


A. Entanglement entropy at the RK point: the Z2 liquid 

One of the most interesting features of the Hamiltonian in 
Eq. (7) is the existence of the exactly solvable RK point {g = 
Erk)- We start by showing that at this point the ground state 
of our model (EPL) indeed has a topological order. 

In order to characterize the topological order, an instructive 
quantity to look at is the entanglement entropy (S') for a bipar¬ 



FIG. 4: (Color online) The plot shows the finite size dependence 
of the topological entanglement entropy 7 as a function of L, the 
perimeter of the subsystem A. The green dashed line corresponds to 
the saturated value of In 2 in the thermodynamic limit. The lattice is 
constmcted on an infinite cylinder assuming periodic boundary con¬ 
dition in one direction. The black dashed line indicates the bipartition 
of the cylinder separating the subsystem A from the rest B. The inset 
shows the linear growth of the entanglement entropy S with L. 

tition of the system into two parts A and B . The entanglement 
entropy of the reduced density matrix pA of subsystem A is 
defined as 5" = —Tr [pA log pa\ where pA is obtained from 
the full density matrix by tracing out all the degrees of free¬ 
dom in the rest {E). For a gapped and topologically ordered 
gapped ground state, S satisfies the “area law” which goes as, 

S{L) =aL-j + O(L-i) -F • • • , (8) 

where a in the leading term is a non-universal coefficient and 
L is the perimeter of the subsystem A. The sub-leading term 
7, also known as the topological entanglement entropy,^’“^° 
is, however, universal bearing the anyonic content of the state 
that reflects the topological order. This is directly related to 
the total quantum dimension (D) of the underlying topological 
field theory as 7 = log D. Since D = 2 for a gapped Z 2 liquid 
(described by the Z 2 gauge theory), the quantity 7 in Eq. (8) 
saturates to log 2 in the thermodynamic limit i.e. L —^ 00. 

The ground state at the RK point can be exactly repre¬ 
sented by tensor networks using the framework of projected 
entangled-pair states (PEPS).^*’^^ We use this construction on 
a cylindrical triangular lattice (Fig. 4) to calculate S' as a func¬ 
tion of L and obtain 7 (Fig. 4) using Eq. (8). The subsystem 
A is constructed by bipartitioning the semi-periodic triangular 
lattice with the dashed line as shown in Fig. 4. The circum¬ 
ference L of the cylinder enters Eq. (8) as the perimeter of the 
subsystem and should be much larger than the maximum cor¬ 
relation length associated with the state.The inset of Fig. 4 
shows the expected linear growth of S with L that is predicted 
by the leading term in Eq. (8). The topological entanglement 
entropy (7) is extracted from the intersection of the function S 
on the y-axis when extrapolated backward and plotted against 
different values of L. The tendency of 7 to saturate at the 
value of log 2 for large L indicates to the fact that the RK 
point for the FPL model on the triangular lattice represents 
the ground state of a topologically ordered Z 2 liquid akin to 
the other dimer models at 1 /6 and 1 /2 fillings.This is one 
of the main results of this work. 
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(a) (6) 





FIG. 5: (Color online) (a) The real space dimer-dimer correlation is 
plotted on a triangular lattice cluster at Vrk = — 3.0 which indicates 
the parallel loop pattern deep into the LN phase, (b) Same is plotted 
for Vrk = 1.0 showing the exponential decay of the dimer-dimer 
correlation function in the liquid phase at the RK point. The refer¬ 
ence link, with respect to which correlations at all the other links are 
measured, is denoted by Ro- Red and blue links stand for positive 
and negative correlation respectively while the width measures the 
correlation strength, (c) Plotted is the same correlator as a function 
of |ri — rj I at different RK potential. The red dashed line represents 
an approximate exponential fit for ^ ~ 1.0. 

However, away from the RK point in the accessible param¬ 
eter space, the ground state is no longer exactly known. We 
adopt ED techniques to infer that the topological liquid is sta¬ 
ble even away from the RK point over a finite window (see 
Fig. 3) up till the system undergoes a quantum phase transi¬ 
tion into the LN phase beyond a critical value of Vrk- In the 
following subsections, we present the ED results containing 
the information about the low-lying spectrum of the model 
and measurements of various correlation functions that reflect 
the salient features of the phase diagram depicted in Fig. 3 
middle line. 


B. Dimer-dimer correlation and the LN order 

At the RK point, the loops are strongly fluctuating and the 
dimer-dimer correlations decay exponentially (as opposed to 
the algebraic decay on the square lattice) like the one in the 
QDM at 1 /6 filling. This indicates the lack of long range 
dimer order, as expected, in the Z 2 liquid phase. To check 
this, we calculate the dimer-dimer correlation function 

- rj\) = (n(r,)n(rj)) - (n(r,))(n(rj)) , (9) 

using ED on different clusters. Here n{ri) is the dimer oc¬ 
cupation number on link of the triangular lattice specified 


by the position vector r^. Each cluster has an extension of 
along the unit vector ai = (1,0) and Ly along the unit vector 
a2 = (5, (see Fig. 2 ). 

Although the considered clusters are small, the numerical 
results shown in Fig. 5(c) suggest that the correlator (as a 
function of r = jr^ — r^l) is indeed exponential at the RK 
point with the correlation length ^ ^ one lattice spacing. This 
is a known result for the case of the QDM on a triangular 
lattice (or equivalently the 1/6 filling).^® As mentioned ear¬ 
lier, at the RK point the model is exactly solvable and all 
the quantum correlations can be calculated classically (using 
some classical numerical techniques^®). The real space real¬ 
ization of the dimer correlation function at the RK point is 
shown in Fig. 5(b) which points to an exponentially decaying 
nature of the correlation function as expected. The red and 
blue bonds correspond to the positive and negative correlation 
respectively with the width of the bond being proportional to 
the absolute value of the correlation function. 

The exponential decay disappears as Vrk is taken to large 
negative values till we get into a phase with long range dimer- 
dimer correlations with C'(r) ~ (up to an offset) deep 

inside this phase signaling the onset of long range dimer order. 
To explore the nature of this long range ordered phase we plot 
the real space dimer-dimer correlation function in Fig. 5(a). 
The plot suggests that the ordered phase is characterized by 
parallel alignments of loops on the triangular lattice which 
does not break any translation symmetry but the three-fold 
rotational symmetry (corresponding to the three-fold quasi¬ 
degeneracy in the ground sate of the spectrum). Since the 
parallel loops do not have an orientation along the direction 
of their alignment, we refer to this as the LN phase. 


C. The phase transition between the Z2 liquid and the LN 

Having characterized the Z 2 liquid phase (by non-trivial en¬ 
tanglement entropy and short range dimer-dimer correlations) 
and the LN phase (by long range dimer-dimer correlation), we 
now focus on the quantum phase transition between the two 
phases. First we study the excitation gap in the liquid with 
ED and observe that the gap closes only at a finite distance 
from the RK point (see Fig. 6). Thus the Z 2 topological liquid 
persists over the whole region from 14k = 1 to Vrk ~ —0.3 
which suggests that alone the kinetic term in Eq. (7) can po¬ 
tentially stabilize the liquid phase even for Vrk = 0. The con¬ 
tinuous vanishing of the gap near 14 k = —0.3 gradually de¬ 
stroys the liquid state driving the system into a charge-ordered 
LN phase as suggested from the behavior of the density- 
density correlation function shown in Fig. 5(c). 

A generic way to locate the critical point within ED is to 
look at the response function of the system as a function of the 
parameters that define the Hamiltonian in Eq. (7). In our case, 
an equivalent quantity can be framed using the second deriva¬ 
tive of the ground state energy with respect to the RK poten¬ 
tial: —S^Eg/SV^Y^. This is plotted in Fig. 7 where the single 
peak in the response is visible approximately at Vrk —0.3. 
This shows that a single transition separates the topological 
liquid and the LN. While it is impossible to rule out the possi- 
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FIG. 6: (Color online) In the left panel, the gap to the first excited 
state tends to vanish at 14 k ~ —0.3 for different symmetric and 
stripped clusters. The ground state always lies in the zero-flux sector 
(0,0) for all of them. In the right panel, the derivative of the gap is 
plotted as a function of 14k to locate the transition point at 14k ~ 
-0.3. 



FIG. 7: (Color online) The ground state susceptibility —S^Eg/5V^Y_, 
in the (0,0) sector indicates a continuous phase transition between 
two phases for different clusters. The transition point is close to 
14k ~ —0.3 as evident from the right panel of Fig. 6. 


bilities of a first order transition from our finite cluster results, 
the smooth and gradual increase of the response along with 
the absence of shoulders suggests that the transition between 
the liquid and the LN may be continuous or weakly first order. 
The full energy spectrum of each of the clusters (not shown) 
suggests that at large negative 14 k, the ground state becomes 
nearly 3-fold degenerate (the quasi-degeneracy attributes to 
finite-size effect) where the three states are formed by super¬ 
position of the three loop patterns allowed by the C 3 symme¬ 
try. We conclude on our numerical results from ED by noting 
that as the system gradually enters the ordered phase crossing 
the transition point, a set of three states in the bottom of the 
spectrum starts separating from the rest. Deep inside the LN 
phase, these three states become nearly degenerate (up to the 
finite size effects) with a finite excitation gap which is much 
higher in magnitude than the liquid gap and scales linearly 
with 14 k- This quasi-degeneracy is exact in the thermody¬ 
namic limit where the C 3 symmetry is spontaneously broken. 


Having established the two phases and the possibilities of 
a continuous phase transition between them, we now explore 
the critical theory for the predicted critical point (at Vrk 
—0.3). We note that such a continuous transition would be 
very interesting in the sense that it describes the destruction 
of a topologically ordered phase. 


IV. CONTINUOUS TRANSITION BETWEEN THE Z2 
LIQUID AND LN 

To construct a theory for the continuous phase transition be¬ 
tween the Z 2 liquid and the LN phases, we now introduce an 
alternative spin representation of the hard-core boson model 
(or the equivalent dimer model). 


A. Spin representation and the gauge theory 

To obtain such a description, we first identify a spin 1/2 
degree of freedom on each site of the kagome lattice by virtue 
of the well known mappings: bl = af, bi = a~ and 
+ l )/2 where an up(down) spin represents pres- 
ence(absence) of a boson at the lattice site. Eq. (4), then, be¬ 
comes 


Hv = V^^(ia£ + /i)'-/r'iVo/16U, (10) 

{O} 

where Nq is the total number of kagome hexagons and h = 
3 — 6 /. The sum of all spin moments in a given hexagon 
is denoted by Sq = where ctq = Clearly 

in the spin description, different fillings of bosons (/) corre¬ 
spond to different integer values of h which essentially plays 
the role of an external magnetic field. Lowest energy con¬ 
figurations of the spin system specified by Eq. (10) satisfy 
the constraint that sum of the moments in every hexagon is 
exactly opposite to h which depends on the filling factor /. 
Thus for / = 1/6,1/3,1/2; h = 2,1,0 and hence poten¬ 
tial term is satisfied if the total magnetization per hexagon at 
these fillings are Sq = — 2 , — 1 , 0 respectively, corresponding 
to having one, two or three up spins (which means presence 
of one, two or three bosons as expected) per hexagon. 


B. Effective spin model 


In terms of the spins, the effective Hamiltonian (in Eq. [5]) 
representing the dynamics within the degenerate ground state 
manifold of 'Hy is given by. 


Heff = ^ Va 

Ol 


-5 n + Trk 

OL 


( 11 ) 


Each term in Eq. (11) involves a product of four spins that 
form a bow-tie, as shown in Eig. 10(a). The projector which 
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FIG. 8: (Color online) The two-vison correlator as a function of | r/ — 
rj| at different values of Vrk for a symmetric cluster of L = 4. 
The red dashed line represents an exponential decay; 
for C ~ 0.7. 

selects out the flippable bow ties is expressed as'^ 

J±1 aGa ^ 

and we have also added the potential term in Eq. (11) to re¬ 
cover the RK Hamiltonian given in Eq. (7). 

When the first term in Eq. (11) dominates, the system 
prefers to align all the spins in the cr® direction and hence the 
boson number per site fluctuates. This is indeed the salient 
feature of a Z 2 liquid phase (see below). On the other hand 
when the second term dominates, the system prefers to choose 
a pattern to order in the direction and we have a long 
range order in the boson density which turns out to be the LN 
phase that we discussed before. The actual magnitude of the 
coupling constants for which the transition between the two 
phases takes place depends on the details of the microscopic 
model. 


C. Ising gauge theory, visons and vison correlator 

Each spin sitting at the site of the kagome lattice is a part 
of two hexagonal plaquettes. We can dehne the Ising vari¬ 
ables e/i = ±1 for each such hexagonal plaquette. Clearly 
the Hamiltonian in Eq. (11) is invariant under the I 2 gauge 
transformation'® 

af ^ , (13) 

where h and h' denotes the two hexagons of which the site i 
is a part of. Such an Ising gauge structure, generated by the 
gauge transformations Qh = exp[i7rS'Q], is an emergent prop¬ 
erty of the low energy subspace*® of the original microscopic 
bosonic Hamiltonian (in Eq. [1]) in the strong coupling limit 
(V ^ t). Indeed the above Hamiltonian represents an Ising 
gauge theory where the plaquette term ) and 

= (crfcrfcrf), shown in Eig. 10(a), measures the Ising 
magnetic flux through each triangle of the kagome lattice. 

The operator erf creates two such magnetic fluxes on the 
two triangles of the kagome lattice of which it is a part of. 


Such Ising flux excitations have been dubbed as visons .One 
can create a single vison excitation'® by applying a product 
of (T* operators along any path C starting from any site of a 
given triangle of the kagome lattice to inhnity; 

I—>^00 

v} = n < ’ (14) 

c 

where i is a kagome lattice site encountered on the path C 
which runs from the corresponding triangular kagome plaque¬ 
tte I, where the vison resides, to spatial 00. Since the path op¬ 
erators commute with each other, it is straightforward to show 
that the two-vison wave function is symmetric under the ex¬ 
change of the visons, or in other words, the visons are bosons 
themselves. With reference to the dimer covering, the above 
operator is nothing but the number of dimer variables encoun¬ 
tered along the path. This immediately implies that the vison- 
vison correlator is given by®*'’®* 

i^j 

=VIJ . ( 15 ) 

c 

(where the path C runs from I to J) translates to 

(vM) = ((- 1 )^-) , ( 16 ) 

where the operator Njj counts the total number of dimers en¬ 
countered on the path C form / to J in any given dimer con¬ 
figuration. It is easy to verify that the explicite expressions of 
the vison operators in terms of the spins given in Eq. (18) are 
independent of the contour C up to an overall sign that can 
be hxed by measuring the correlator with respect to a hxed 
reference dimer conhguration. 

In the Z 2 liquid phase, where J^a = = +1’ 

son excitations have a hnite gap and the ground state does not 
contain free visons. On the other hand, the phase character¬ 
ized by (tr^) ^ 0, which we shall show is the LN phase, is a 
vison condensate. Thus we expect that the transition between 
the Z 2 liquid and the LN is described by the closing of the vi¬ 
son gap leading to the vison condensation. Hence, the ground 
state expectation value of the vison correlator should also de¬ 
cay exponentially in the liquid phase with a length scale pro¬ 
portional to the inverse of the vison gap while it should have 
asymptotically reached a constant value in the LN phase. 

In Eig. 8 we show the two-vison correlation function at 
different values of Vrk- The data at Vrk = 1.0 hts well 
with the exponential curve for ^ ~ 0.7. As Vik 

is decreased further, the system gradually enters the ordered 
phase and the two-vison correlator becomes asymptotically 
constant. This behavior is expected, as we shall show below. 
Note that in Eig. 8, close to Vrk = —0.3 the value of vjj 
changes by an order of magnitude much in the similar way as 
the density-density correlation does in Eig. 5. 

D. Lattice description for the visons 

Our numerical results indicate that the phase transition be¬ 
tween the Z 2 liquid and the LN phase, driven by the conden¬ 
sation of vison excitations, is possibly continuous. We shall 
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now derive an effective critical theory for such a continuous 
transition. This would then compliment our numerical under¬ 
standing of the phase diagram of the microscopic model. 

In spirit of the universality of continuous phase transitions, 
to this end, we perform a series of mappings to isolate the 
vison degrees of freedom which we use to describe the critical 
theory for the transition. The effective Hamiltonian in Eq. (11) 
can be obtained from 

+ -^v E (^O + 

a O 

+ UeS Vg (17) 

a 

in the limit of Peff/Kff, Weff/Ks 0 Kff > 0) where 
the leading term is obtained in the second order perturbation 
theory with g ~ Pgjf/14ff and Ueff ^rk- We immediately 
note that the last two terms in Eq. (17) commute with each 
other as expected and hence, in the regime |ueff | 3> Peff (and 
VeB being the largest energy scale), the above model is ren¬ 
dered classical and the Ues term chooses an appropriate order 
that is consistent with the filling (i.e. allowed by 14ff). Eor 
Meff < 0, this favors the order as depicted in Eig. 10(c). On 
the other hand, when |ueff| ^ geS, this order is expected to 
melt due to the quantum fluctuations in the spins. In particu¬ 
lar, a state which allows flipping of spins in closed loops while 
maintaining the Ailing constraint becomes favorable. This is 
indeed consistent with the picture of the Z 2 liquid of the mi¬ 
croscopic model which is exact at the RK point. 

To arrive at the critical theory describing the vison conden¬ 
sation, we exploit the gauge structure, described in the last 
sub-section, of the effective theory to isolate the vison degrees 
of freedom. At this point, we introduce the honeycomb lattice 
which is the medial lattice of the kagome lattice as shown in 
the left most panel in Eig. 9. On the sites of this medial lattice, 
we define the Ising variables = ±1 and on its links we de¬ 
fine the Ising gauge fields = ±1. The mapping from the a 
to the V and p variables is defined as follows, 

= < = vfphvy ( 18 ) 

A 


Clearly, with the insight of Eq. (14), in the above mapping the 
ufs are nothing but the visons (with conjugate momenta given 
by vf). Since they carry Ising magnetic charge, they cou¬ 
ple to the dual Ising gauge potential, pjj and transform under 
a projective representation (that forms a projective symmetry 
group) of the symmetries of the underlying spin Hamiltonian. 
We note that the product of the dual gauge fields around the 
hexagonal plaquette is given by 

l[ph = l[<=^^^ (19) 

o o 

where the first product is over each honeycomb hexagon and 
the second one is on each kagome hexagon. 

Though not directly relevant to this work, we would like 
to point out that the presence of the dual magnetic flux, is 
nothing but the Ising electric charges (spinons): 


n'>?./=ii'’:= 

o o 



FIG. 9: (Color online) The medial honeycomb lattice is shown along 
with the original kagome lattice. Removal of the blue-colored bo¬ 
son creates two spinon excitations sitting at the centers of the ad¬ 
jacent shaded hexagons (the “defect hexagons”) which violate the 
constraint of having two bosons. By virtue of the single particle hop¬ 
ping term, these two defects can separate further apart, however their 
creation is energetically suppressed in the strong coupling limit. 


which sit at the center of the hexagonal plaquettes. In terms of 
the original bosons, these represent hexagons where the con¬ 
straint of having two bosons is violated (see Eig. 9). Evidently, 
such “defect hexagons”, are created in pairs and are energeti¬ 
cally very costly in the strong coupling limit. It is clear from 
the above equation that the visons see the spinons as a source 
of TT-flux and hence the dual gauge potential naturally cap¬ 
tures the mutual semionic statistics between the visons and the 
spinons.'® Now note that at 1/3 filling, we have two bosons 
per hexagon. Hence, we must have 


= = ^ ( 21 ) 
o o 

This just means that there are no spinons at low energies be¬ 
cause they are too costly. Hence, on circling the plaquettes of 
the medial honeycomb lattice, the “flux” seen by the u^-spins 
is zero. This type of Ising gauge theory (IGT) are called even 
IGT as opposed to the odd IGT that arises in case of 1 /6 or 
1 /2 Ailing when the u^-spins sees a tt-Aux in each honeycomb 
plaquette. 

With the mapping in Eq. (18) the Hamiltonian in Eq. (17) 
becomes (we put Ues = 0 for the moment) 

Tfdual = ~geS ^ ^ Vj -f Vgff ^ ^ | — ^ ' ^iPlJ^J + ^ | 

I O \ (/,J)GO / 

( 22 ) 


expanding which we get. 


"Hdual = -5eff E E 

I (I,J) 



E ^iPn'vl + E ^iPhv>i'Prj'v}' 

{{II')) i,J,i',J'eo 


(23) 


up to a constant. Until now we have ignored the potential term 
in the RK Hamiltonian. This term, in the spin language, has 
the form 

Ct 


- 1 , 


( 20 ) 


(24) 
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where a refers to all bow-ties and has the typical form, for the 
bow-tie in Fig. 10(a), of 


/O _ r ^ 2^2 ^ 2^2 ^ 2^2 i ^ 2^2 

^OL — g [ <^ 1^2 *^ 3^4 *^ 2^3 ' ^1^3 

^Z] 

'^1^2 ‘^3^41 


^2^4 


(25) 


This, under the mapping to the vj variables, augments the sec¬ 
ond neighbor as well as the four-spin terms along with provid¬ 
ing the interactions between four spins, two each on adjacent 
hexagons. As remarked earlier, the fields transform under a 
projective symmetry group (PSG) and in general a third neigh¬ 
bor Ising term for the vf fields would also be allowed by the 
PSG of an even-IGT. These term, would arise, for instance, 
on integrating out the four spin interactions which would also 
renormalize the nearest and the second neighbor interactions. 
In addition, a Maxwell term for the pjj fields of the form 
rio PiJ allowed and they renormalize the energy of 

the spinons, but due to the constraint in Eq. (19), such terms 
are trivial and hence left out. Now, Because of the constraints 
in Eq. (19) {i.e. no spinons), it is possible to choose a gauge 
where 


p]j = +i, yi,j. ( 26 ) 

Hence the general form of the model is given by 

^dual = JlY. + + E ■ 

n.n n.n.n n.n.n.n J 

(27) 

Thus, the minimal gauge theory for the EPL model is dual to 
the ferromagnetic transverse field Ising model on the honey¬ 
comb lattice with first, second and third neighbor Ising inter¬ 
actions. 

In the large T limit we can neglect the J terms and the v- 
spins are polarized in the x direction with the finite vison gap 
^ 2r. The underlying spins are thus fluctuating allowing 
for the boson number per site to fluctuate. This paramagnetic 
phase is nothing but the Z 2 liquid in the dual description. 

On increasing the Ising couplings (^ 1 , 2 , 3 ), the visons gain 
dispersion and if the minima of the dispersion touches zero, 
they can condense leading to (uf) ^ 0. The nature of the 
ordering depends on the relative signs and magnitude of the 
Ising couplings and the phase diagram in the limit of large 
> 7 i, 2 , 3 /r is shown in Eig. 10(d). Setting J 2 /J 1 = G and 
J3/J1 = t 2 and r / Ji = f with Ji = 1, in Eq. (27), the soft 
vison modes^^’® can be obtained from the Eourier transform 
of the Ising terms. This is given by 


Tfdual — E J{k)^ly>k , (28) 

k 


where Jik) is the Eourier transform of the adjacency matrix 
of the Ising terms in Eq. (27): 


J{k) = 


f 2tiS 

Vt* + t 2 V* 


7 + hv\ 
2 tiS J 


(29) 



FIG. 10: (Color online) (a) A flippable bow tie. (b) The first Brillouin 
zone of the honeycomb lattice showing M points, (c) Loop ordering 
at Ml point. A loop segment is placed on the link of the triangular 
lattice whenever it crosses a honeycomb bond joining antiparallel 
spins which represent the vison degrees of freedom, (d) The yellow 
region qualitatively shows the allowed ranges of ti and t 2 for which 
the minima of J{k) lie in the vicinity of the M points within an 
energy window of 10“^. The vertical phase boundary indicates to 
the case when the minima are exactly on the M points corresponding 
to the LN phase for Vrk <C 0 . 


and vpit = and The momen¬ 

tum dependence of the parameters goes as 


(5 = cos fci-F cosfc 2 -F cos (fci-F fc 2 ) 

7 = 1 + 

T] = 2 cos (/ci-F ^ 2 ) + (30) 

where ki = k ■ st with si = (•\/3/2, —3/2) and 5*2 = 
(•\/3/2, 3/2) being the basis vectors of the honeycomb lat¬ 
tice shown in Eig. 10(c). Eor an extended range of positive 
ti and t 2 , as highlighted by the yellow region in Eig. 10(d), 
the minima of the energy dispersion J{k) occur at the three 
inequivalent M points; Mi, M 2 and M 3 in the Brillouin 
zone [see Eig. 10(b)] whose coordinates are (tt/v/S, - 7 r/ 3 ), 
(tt/v/S, 7 r/ 3 ) and (0, 27r/3) respectively. This extended re¬ 
gion of the phase diagram, where the ordering occurs at the M 
points of the Brillouin zone of the medial honeycomb lattice, 
yields to a u| ordering pattern which is shown in Eig. 10(c). 
Translating back (using Eq. (18) and Ui = {1 + af)/2), this 
gives rise to the ordering pattern for 14 k ‘C 0. There are three 
such patterns (for the three M points) which corresponds to 
the three LN phases related by the C 3 symmetry as obtained 
from our previous numerical calculations. 

While it is not possible to calculate the values of the Ising 
couplings in terms of the couplings of the underlying RK 
Hamiltonian in Eq. (7), we observe that, the effective model 
allows phases that are observed in the microscopic model to¬ 
gether with the possibilities of a direct continuous quantum 
phase transition from the Z 2 liquid to the LN phase. Since 
the symmetries of the actual QDM and the effective gauge 
theories are identical, we expect that the transitions, which is 
attributed to the condensation of the visons, in both the mod¬ 
els belong to the same universality class. We now attempt to 
And the structure of the critical theory which can predict the 
nature of this transition. 








10 


E. Critical theory for the transition 

Approaching the transition point from the liquid side ( F > 
J in Eq. [27]), we can write down the critical theory in terms 
of the vison modes that goes soft at the transition. These soft 
modes can be written as®' 


3 

= (31) 

i=i 

where are the amplitudes of the three soft modes occurring 
at the three M points of the Brillouin zone. To construct the 
Landau-Ginzburg action in terms of the soft modes we need to 
figure out the transformation of il^j (j = 1, 2, 3) among them¬ 
selves under various symmetries of the Hamiltonian. These 
are: ( 1 ) Ti: lattice translation along si, ( 2 ) T 2 : lattice trans¬ 
lation along S 2 , (3) I: bond inversion or parity, (4) Cg: rota¬ 
tion of 7 r /3 about the center of a plaquette, and (5) global Z 2 
symmetry under which —v^. On a point of specific co¬ 

ordinates {x,y}, the actions of the above symmetries are the 
following: 

Ti : {x,y]a,b} ^ {x + l,y\a,b} 

T 2 : {x,y]a,b} ^ {x,y + l\a,b} 

I : {x,y]a,b} ^ {-x,-y\b,a} 
\{x,y\a}^{x-y + l,x]h}, 

^6 : r , T r T ( 42 ) 

{x,y\b\ {x - y,x;a|. 


The transformation matrices of the three critical modes cor¬ 
responding to different lattice symmetry operations and the 
global Z 2 are as follows: 



■-1 

0 

0 



'l 

0 

0 ■ 

Ti-Ti = 

0 

1 0 


to 

II 

0 

-1 

0 


0 

0 - 

1 



0 

0 

- 

1 


'1 0 

0 ■ 




0 

0 

-1 


TZj = 

0 1 

0 


5 


1 

0 

0 



0 0 

-1 




0 

1 

0 



■-1 

0 

0 ■ 






T^I.2 = 

0 

-1 

0 







0 

0 

- 

1 







These five matrices generate a 24 elements finite subgroup of 
0(3) which is isomorphic to C 2 ® 2 ) 4 .®^ Respecting all the 
above projective symmetry transformations, the most general 
Landau Ginzburg (LG) functional in (2 -f 1) dimensional Eu¬ 
clidian space-time assumes the following form. 


5 = 


d^vdT C , 


(34) 


where the Lagrangian density (up to 6 th order) is given by 

£ = • 's/'ijj + dr'if ■ drip + rip ■ ip + u{tp ■ ipY 

+ v{ip ■ ipf + aiipl + V '2 + ipt) + (35) 


with Ip = ('i/ji, '!/; 2 , If a = 6 = 0, Eq. (35) represents 

the usual soft spin 0(3) action (or the N = 3 linear a model). 
The a term introduces cubic anisotropy in the system. 

Evidently, for a < 0, the ordering occurs at one of the soft 
modes preferentially, i.e., the functional is minimized when 
one of the components (among ipi, ip 2 , and ip'^) takes a finite 
value while other two remain zero. Three such possibilities 
give rise to three symmetry oriented loop orderings. Eor ex¬ 
ample, the order at Mi, which corresponds to ^ 0 and 
^^2 = tp 3 = 0 ^ can be read from the structure of the full eigen¬ 
vector vi = (1,1)^ of J{k) in Eq. (29) set to the momentum 
Ml = (tt/TS, — 7 r/ 3 ). The resultant u^-spin configuration is 
shown in Eig. 10(c) and the loop covering of the dimers can be 
obtained using Eq. (18) which is equivalent to replacing the 
antiferromagnetic bonds of the honeycomb lattice by a loop 
segment on the underlying triangular lattice. This gives back 
the loop ordering which does not break any translation sym¬ 
metry (because Mi and —Mi are identical and related by the 
momentum space lattice vectors) but spontaneously breaks the 
rotational symmetry. Similarly the ordering at other M points 
are related to the present one by the spontaneously broken C 3 
symmetry. This way the critical theory captures the ordering 
patterns obtained in the numerical calculations of the micro¬ 
scopic model. 

Eor 6 = 0, leading order e = 4 — d expansion suggests 
that the cubic anisotropy is irrelevant and the critical point is 
of 0(3) Wilson-Eisher type®^*®® (however, higher order ex¬ 
pansion suggests that it may belong to the cubic critical point 
with critical exponents very close to the 0(3) class®®). The 6 - 
th order anisotropy term (denoted by 6 7 ^ 0 ) is also irrelevant 
at such a point along with the 0(3) invariant 6 th order term, 
V in Eq. [35]. These considerations suggest that the phase 
transition between the Z 2 liquid and the LN phase, as seen in 
the microscopic model, may belong to the 0(3) universality 
class and the anisotropy terms are dangerously irrelevant at 
this critical point. 

It is important to note that the critical theory is not written in 
terms of the order parameter of the LN phase, as the conven¬ 
tional theory of phase transition would suggest.®^’®® Instead, 
it is naturally written in terms of vison fields and the order 
parameter is bilinear in terms of such vison fields. Hence we 
should expect large anomalous dimensions in the scaling di¬ 
mension of the LN order parameter.®®’’® Such large anoma¬ 
lous dimensions are characteristics to these types of uncon¬ 
ventional phase transitions. 


V. SUMMARY AND OUTLOOK 


In summary, we have studied the strong coupling limit of 
the extended Hubbard model for hard-core bosons with short 
range repulsive interactions. Eocusing on the particular frac¬ 
tional filling of 1 /3, we show that the low energy physics of 
the bosonic model is described by a quantum EPL model on 
the triangular lattice. Using a combination of different numer¬ 
ical techniques we analyze the ground state phase diagram of 
this quantum EPL model. Our numerical calculations conclu¬ 
sively establish the presence of a topologically ordered Z 2 hq- 
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uid phase over an extended parameter regime of the effective 
low energy Hamiltonian. On tuning appropriate parameters of 
the effective model, we find indication for a continuous phase 
transition from the topological phase to a phase that breaks 
the C 3 rotational symmetry of the triangular lattice. Using a 
mapping to an Ising gauge theory and PSG based arguments, 
we construct the effective field theory for a generic continuous 
transition between the two phases, which we argue, belongs to 
the 0(3) universality class. We then show that such a transi¬ 
tion is naturally related to the condensation of the (bosonic) 
magnetic charges of the Ising gauge theory, the so called vi- 
sons. Hence, contrary to the conventional theory of continu¬ 
ous transition,®^ the present critical theory is written in terms 
of the soft vison modes and not the LN order parameter which 
is a bilinear in terms of the vison fields. 

The present calculations show that the model of strongly 
interacting hard-core bosons can harbor rich and interesting 
phase diagrams including conventionally as well as topolog¬ 


ically ordered states at different fractional fillings. It is in¬ 
teresting to note that considering the interacting particles to 
be fermions would induce a non-trivial statistics in the dimer 
problem. Interestingly, such kind of fermionic dimer mod¬ 
els can arise in context of describing the metallic state of the 
hole-doped cuprates at low hole densities.^' Whether a sta¬ 
ble liquid phase can still be realized in such fermionic models 
requires more understanding regarding the underlying gauge 
theory and constitutes an interesting direction for further stud¬ 
ies in future. 
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